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Abstract 

We show that the Painleve VI equation has an equivalent form of the non- autonomous 
Zhukovsky-Volterra gyrostat. This system is a generalization of the Euler top in C'^ and 
include the additional constant gyrostat momentum. The quantization of its autonomous 
version is achieved by the reflection equation. The corresponding quadratic algebra gener- 
alizes the Sklyanin algebra. As by product we define integrable XYZ spin chain on a finite 
lattice with new boundary conditions. 
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1 Introduction 

In this paper we discuss a few issues related to isomonodromy problems on elliptic curves, 
integrable systems with the spectral parameter on the same curves and the XYZ spin-chain on 
a finite lattice. Our main object is the Painleve VI equation (PVI). It is a second order ODE 
depending on four free parameters (a, /?, 7, S) 

d?X _l f I 1 1 \(dxV fl 1 I \dX 

IF ~ 2\X^ X - l^ X - t) \ dt J ~ t-l^ X -t) 'dt^ 

X(X-l)(X-t)f t-l t(t-l)\ , , 

PVI was discovered by B. Gambler ^HI in 1906. He accomplished the Painleve classification 
program of the second order differential equations whose solutions have not movable critical 
points. This equation and its degenerations PV — PI have a lot of applications in Theoretical 
and Mathematical Physics, (see, for example j41j). 

We prove here that PVI can be write down in a very simple form as ODE with a quadratic 
non-linearity. It is a non-autonomous version of the SL(2, C) Zhukovsky-Volterra gyrostat (ZVG) 
|42| I39j . The ZVG generalizes the standard Euler top in the space by adding an external 
constant rotator momentum. The ZVG equation describes the evolution of the momentum 
vector S = {Si, 82,33) lying on a SL(2,C) coadjoint orbit. We consider Non-Autonomous 
Zhukovsky-Volterra gyrostat (NAZVG) 

drS = Sx (J(t) • 5) + 5 X z7' (1.2) 

where J(r) • S = {JiSi, J2S2, J3S3). Three additional constants u' = (i/[, 1/3, 1^3) form the 
gyrostat momentum vector, and the vector J = { Ja(T)}, (a = 1, 2, 3) is the inverse inertia vector 
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depending on the "time" r in the following way. Let S,- = C/(Z + rZ) be the elliptic curve, 
and p{x,t) is the Weierstrass function. Then Ja(T) = p{uJa,T), where coa are the half-periods 
(l' i' ^^^)- "^^^ constants Ua along with the value of the Casimir function S"^ = (^q)^ are 
expressed through the four constants (a, /3, 7, 6) of the PVI. If the gyrostat momentum vanishes, 
z?' = the ZVG is simplified to the non- autonomous Euler Top (NAET) equation 

drS = S X {J{t) ■ S) . (1.3) 

To establish the connection between the PVI eq. and (|1.2|) we start with the elliptic 

form of (mU) Hi HI 

(fu ^ 

-j^ = -'^^lp'{u + U;a,T) , {uJa = iO,UJa)) ■ (1-4) 
a=0 

The interrelation between two sets of the constants Va and is explained in Sections 5.3 and 
5.4. We call this equation EPVI to distinguish it from By replacing r on the external 

time t we come to the BCi Calogero-Inozemtsev system (CI) 

= -^^lp'{u + ^a,T). (1.5) 
a=0 

This equation is the simplest case of the integrable BCat CI hierarchy When the constants 
are equal, ()1.5|) describes the two-body elliptic Calogero-Moser (CM) system in the center of the 
mass frame 

^ = -z.V(2n,T). (1.6) 
In this case EPVI (|1.4)) assumes the form 

^ = -u^dup{2u,T). (1.7) 

The genuine interrelations between integrable and isomonodromy hierarchies arise on the 
level of the corresponding linear systems. The linear equations of the isomonodromy problem 
look like a quantization (the Whitham quantization) of the linear problem for the integrable 
hierarchy |2nj . While for the integrable systems the Lax matrices are sections of the Higgs 
bundles (one forms on the spectral curve), they become the holomorphic components of the flat 
connections for the monodromy preserving equations. Recently, the Lax representation for (|1.5|) 
was proposed in Ref. |42j. It allows us to construct the linear system for EPVI (|1.4|1 . 

In our previous work we proposed a transformation of the linear system for the N-body 
CM system to the linear system for the integrable SL(A^, C) elliptic Euler- Arnold top 30 . It is 
the so-called Hecke correspondence of the Higgs bundles. It is accomplished by a singular gauge 
transform (the modification) of the Lax equations. This action on the dynamical variables is a 
symplectomorphisms. In the simplest case it provides a change of variables from the two particle 
CM JUni) to the sl(2, C) autonomous Euler Top (ET) 

dtS = Sx{J-S), {J ■ S) = (JaSa) . (1.8) 

It opens a way to define the Lax matrix L^'^{z) of the ET from L^^^ . The both models contain 
a free constant which is the coupling constant for CM and the value of the Casimir function of 
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the SL(2,C) coadjoint orbit for ET. It should be mentioned that the analogous transformation 
was used in ^HHHIISHI for other purposes. 

In the similar way we prove that the autonomous SL(2, C) ZVG 

dtS = S X JS + S xiJ' . (1.9) 

is derived from the BCi CI eq. (|1.5|) via the same modification. In this way we obtain the 
Lax matrix L^^^{z) from the Lax matrix of the CI system. The change of variables is given 
explicitly. For SO (3) ZVG the Lax pair with the rational spectral parameter was constructed 
in Ref. [H] and discussed in Ref. In our case the Lax pair depends on the elliptic spectral 
parameter. 

In fact, the modification can be applied to the isomonodromy problem |^|^. It acts on 
connections and, in particular, transforms (|1.7() to NAET H1.3() and the generic EPVI (|1.4|) to 
the NAZVG (|1.2j) . We present the explicit dependence S{u,dru)- In this way we establish the 
equivalence between PVI (dHJ and NAZVG ((T3|) . 

There exists another way to define the Lax matrix L^^'^{z) and thereby to derive ZVG 
H1.9|) that will be used in this paper. The starting point is a special Elliptic Garnier-Gaudin 
system (EGG). EGG is an example of the Hitchin systems |22]. It is derived from the rank two 
quasi-parabolic Higgs bundle [SJ over S,- of degree one. We assume that the Higgs field has 
simple poles at the half-periods uja, {a = 0, ... ,3). The invariant part of the Higgs field with 
respect to the involution z — > —z leads to L^^'^{z). 

Finally, we can derive ZVG starting with the quantum reflection equation. Consider flrst 
the quantization of SL(2, C) ET. It can be performed by the quantum exchange relations with 
the Baxter i?-matrix ^3': 

R{z, w)Lf^{z)Lf^{w) = L^^{w)Lf^{z)R{z, w) , (1.10) 

where L^'^{z) is the quantum Lax matrix of the elliptic top. This equation is equivalent to the 
Sklyanin algebra j34j . We show that the quantum Lax matrix for ZVG satisfies the reflection 
equation introduced in jS^j: 

R-{z, w)Lf^^{z)R+{z, w)L^^^{w) = L^^^{w)R+{z, w)Lf^^{z)R-{z, w) . (1.11) 

The corresponding quadratic algebra generalizes the Sklyanin algebra. In the classical limit 
(|1.11() yields two Poisson structures for ZVG: 

{Li{z),L2{w)}2 = ^[L,{z)L2{w),r-{z,w)] - h^{z)r+ {z,w)L2iw) + h2{w)r+ L^{z) , (1.12) 

{Li{z),L2iw)}i = ^[Li{z)+L2H,r-{z,w)] - ^[Li{z) - L2iw) , r+ {z , w)] 

which are compatible as in the case of the Sklyanin algebra J7]. The first type of brackets 
generalizes the classical Sklyanin algebra while the second is just Poisson-Lie brackets. The 
coefficient ^ in (|1.12|) comes from the statement that these brackets are derived from the standard 
brackets 

{Li{z),L2{w)}2 = [Li{z)L2{w),r-{z,w)] , (1.13) 

{Li{z),L2{w)}i = [Li{z) + L2{w),r-{z,w)] 

by the Poisson reduction procedure for the constraints L{z) + L^^ {— z) det L{—z) = and 
L{z)+L{—z) = for the quadratic and linear brackets correspondingly. This procedure however 
will not be discussed here. 
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In jnS] the reflection equation was used to construct an integrable version of the XYZ spin- 
chain on a finite lattice. Following this recipe we obtain the XYZ spin-chain with the quantum 
ZVG on the boundaries. Three additional constants here combine into the vector of magnetic 
field. The classical Hamiltonian is presented in Proposition 8.3. As far as we know the obtained 
model was not discussed earlier. 
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2 Isomonodromic deformations and Elliptic Calogero-Moser 
System 

We consider differential equations related to the A^-body integrable elliptic Calogero-Moser sys- 
tem with spin (CM) [HI I24| 111! I4()j . They are defined as monodromy preserving equations of 
some linear system on an elliptic curve and generalize to A^ dependent variables 20 . It is 
a Hamiltonian non-autonomous system that describes dynamics of A^ particles with internal de- 
grees of freedom (spin) in a time-depending potential. We call this system Non-autonomous 
Calogero-Moser system (NACM). In section 2.4 we consider the interrelations between au- 
tonomous equations, corresponding to integrable hierarchies and non-autonomous (monodromy 
preserving) equations. 

2.1 Phase space of NACM 

The phase space of NACM system is the same as of CM. Let = C/zf^, 'zf^ = Z (B tZ, 
(9m r > 0) be the elliptic curve. The coordinates of the particles lie in S^: 

U = {Ul, . . . ,Un) , Uj G T.r 

(2) 

with the constraint on the center of mass 

J2 Uj = 0. Let Z)- ' X Wn be the semi-direct product 
of the two-dimensional lattice group and the permutation Wn ■ The coordinate part of the phase 
space is the quotient 

A = (C^/(Z(2) xt^^))/S.. (2.1) 
The last quotient respects the constraint on the center of mass. Let 

V = {vi,. . . ,vn) , vjeC, y^^vj = 0. 

be the momentum vector dual to u : {vj,Uk} = 6jk- The pair (v,u) describes the "spinless" 
part of the phase space. 

The additional phase variables, describing the internal degrees of freedom of the particles, 
are the matrix elements of the A^-order matrix p. More exactly, we consider p as an element 
of the Lie coalgebra sl(A^, C)*. The linear (Lie-Poisson) brackets on sl(A^, C)* for the matrix 
elements have the form 

{Pjkj Pmn} — Pjn^km Pmk^jn • (^•^) 
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Let O be a coadjoint orbit 



O = {p G sl{N, C)* I p = Adlp", h £ SL{N, C) , p" G Sj*} , 



(2.3) 



where Sj is the Cartan subalgebra of sl(A^, C). The phase space TlP^^ = T*A x O is a symplectic 
manifold with the symplectic form 



uj= {dvA du) - {p^dhh^^dhh-^) , 



(2.4) 



where the brackets stand for the trace. The form lo is invariant with respect to the action of the 
diagonal subgroup D = exp^: h — > hhi, hi G D. Therefore, we can go further and pass to the 
symplectic quotient 

d = 0//D. (2.5) 

It implies the following constraints: 

i) the moment constraint pjj = 0, 

ii) the gauge fixing, for example, as Pjj+i = Pj+ij- 

Example. Let p'' = z^diag(A^ — 1, —1, . . . , —1). Then dimO = 2N — 2. It is the most degenerate 
non-trivial orbit. It leads to the spinless model, since in this case dimO = 0. We should 
represent p*^ in the special form that takes into account the moment constraint (i): 



/ 1 1 
1 1 



P 



iC 



1 1 



1 \ 

1 



1 

0/ 



(2.6) 



For N = 2 these orbits are generic. 

In this way we come to the phase space of the CM 



■rpCM 



{T*{A) X O}, 



(2.7) 



Note that 



dim(7^^g^) = 2iV - 2 + dim O - 2 dim{D) = dim O . 



2.2 Equations of motion and Painleve VI 

The CM Hamiltonian has the form 

1 ^ 

H'^^'^P'^ = - ^ ^1 - ^PjkPkjE2{uj - Uk;r) , (2.9) 
j=i j>k 

where E2{x;t) = p(x;r) + 2r/i(T) is the second Eisenstein function ()A.4|) and r plays the role 
of time. ^ For the orbit, corresponding to ()2.6|) . the spinless Hamiltonian is 

1 ^ 

^CM ^ i ^ „2 _ ^2 ^ (2.10) 

_J=1 j>k 

^In what follows we replace the Weierstrass function p(x; r) used in Introduction by the Eisenstein function 
E2{x; t). It does not affect the equations of motion. 
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For general non- autonomous Hamiltonian systems it is convenient to work with the extended 
phase space by including the time. Here we deal with the extended space 7^^^'* = {TZ^^,T), 
where T = {r € C | 9mr > 0}. Equip it with the degenerate two- form 



uj^^t = uj- -d/7^^^''*P^"(v, u, r) A dr , 



where k > is the so-called classical level. Note that w^^* is invariant with respect the modular 
transformations PSL2(Z) of T It means that w*^^* can be restricted on the moduli space 
M = T/PSL2(Z). The vector field 

j,k,l \ ^ ^ rnn / 

annihilate lo^^^ and define the equations of motion of NACM system 

d/(v,u,p,T) 

= Vr/(v, u, p, r) . 

In particular, 

ndrUj = Vj , (2-11) 

KdrVn = - y^^PjkPkjdu„E2{uj - M„; t) , (2.12) 

Ka,p = 2[Ju(T) •p,p], (2.13) 
where the operator Ju • p is defined by the diagonal action 

Ju(t)-P : Pjk E2{uj - Uk;T)pjk , (2.14) 

i.e. each matrix element pjk is multiplied on E2{uj — u^^t). For = 2 we put ui = —U2 = u, 
vi = —V2 = V = ndrU and come to the two body NACM model 

dlu = -v'^duE2{2u) . (2.15) 

It coincides with (|1.7|) since dup{u) = duE2{u). 

2.3 Lax representation 

The goal of this subsection is the Lax representation of (|2.11j) - (|2.13|1 j2()j . 
2.3.1 Deformation of elliptic curve 

Let = {(x,y) E M|x,y € M/Z} be a torus. Complex structure on is defined by the 
complex coordinate 

= {z = x + Toy , To > 0} , , ~ C/(Z + tqZ) , 

or by the operator annihilating the one form dz. 

Consider the deformation of the complex structure that preserves the point z = 0. Let 
x{z, z) be the characteristic function of a neighborhood of z = 0. For two neighborhoods W D U 
of z = define the smooth function 
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Consider a new complex coordinate w 



w = z {z - z){l - x{z,z)) , (p = To-ro). (2.17) 

P 

The new coordinate defines the deformed elliptic curve S,- = {w = x + ry}. The partial 
derivatives with respect to the new coordinates {w,!])) are 



" TI{^-z + ^Ji^z) 



(2.18) 



where 



a = l + - — —{l-x{z,z)), 13 = -- — —{l-x{z,z))' 
P P 



and jj: is the Beltrami differential 

(3 T-Tof I - x{z,z) - {z - z)dx{z,z) 



P- 



a p yi + l-pi(l-x{z,z)-{z-z)dx{z,z))^ 



Note, that 

P 



0, zeU 



p 

The partial derivatives form the basis dual to the basis of the one-forms {dw, dw). The common 
prefactors in the r.h.s. in H2.18() is irrelevant in our construction. For this reason we replace 
the partial derivatives by the vector fields that only annihilate the corrseponding one-forms 
{dw,dw). We preserve the same notations for them: 



(2-19) 

dn, = dz + pdz , 



The pair {t — tq^t — tq) plays the role of local complex coordinates in the moduli space A4 
of elliptic curves with a marked point in a neighborhood of tq. 

In what follows we replace the complex coordinates {w, w) by the independent coordinates 

w = z---^{z-z){l-x{z,z)), 
w = z , 

and the corresponding independent vector fields {dw,du}), annihilating dw and dw 



dw = dz + pdz . 

Note that dyj = d^, while is independent on fi now. We pass from {w,w) to the chiral 
coordinates {w, w = z) because only the holomorphic dependence on the moduli M is essential 
in our construction. The holomorphic coordinate (r— tq) plays the role of time in the Hamiltonian 
systems we consider here (see (|1.2|) and (|1.4j) ). 
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2.3.2 Flat bundles of degree zero 

Let be a flat vector bundle of rank N and degree over the deformed elliptic curve $],-. 
Consider the connections 

(W, W,T) , 



L^°Hw,w,t), L^^\w,w,t) G C7~Map(S^,sl(iV,C)). 
The flatness of the bundle V^y means 

O^L^o) _ ^g^im ^ [^(0)^^(0)] ^ o_ (2.21) 

By means of the gauge transformations f{w,w) G C°°Map (S,- — > GL(A^, C)) 

the connections of generic bundles of degree zero can be put in the following form: 

1. Z(°) = 0. Then from the flatness I^^Ti} we have 

dn,L^°Hw,w) = 0. 

2. The connection of generic bundles of deg(V^) = have the following quasi-periodicity: 

(u- + 1) = L W [w] , {w + t) = e(u)L(0) (u')e(-u) , 

where the diagonal elements of e(u) = diag(exp(27riui), . . . , exp(27rzMjv)) define the moduli of 
holomorphic bundles. We identify u with the coordinates of particles. In fact, uj, j = 1, . . . , N 
belong to the dual to S,- elliptic curve {the Jacobian), isomorphic to $],-. 

3. We assume that L^'^^ has a simple pole aX w = and 

Res\y,=oL^^Hw) = p. 
The conditions 1, 2, 3 fix L^") up to a diagonal matrix P 

L(o) = P + X, (2.22) 

P = diag(ui, . . .,vn) , 

X = {Xjk} , {j / A;) , Xjk = Pjk<t>{uj -Uk,w) . 

The function <f> is determined by ()A.8|) . The quasi-periodicity of L^^^ is provided by HA.17|) . The 
free parameters v = {vi, . . . , vn) of P can be identified with the momenta. 

The fiatness of the bundle upon the gauge transform amounts the consistency of the system 

a. d^^ = o. 

To come to the monodromy preserving equation, we assume that the Baker-Akhiezer vector ^ 
satisfies an additional equation. Let 3^ be a monodromy matrix of the system ()2.23|1 correspond- 
ing to homotopically non-trivial cycles ^ — > ^'3^. The equation 

Hi. {ndr + mW(u'))^' = (2.24) 

means that dry = 0, and thereby the monodromy is independent on the complex structure of 
r^. The consistency of i. and iii. is the monodromy preserving equation 

a,LW - a^MW - -[lW,mW] =0. (2.25) 
In contrast with the standard Lax equation it has additional term dwM^'^\ 
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Proposition 2.1 The equation \2. g,5)) is equivalent to the monodromy preserving equations 
IF771) . IE71) . irm for lOl) and 

M(°) = {r,fc}, (jy^fc), diagM(0) = 

^jfe = Pjkfiuj - Uk,w) , /(ti, w) = du(t){u, w) . 

Proof is based on the Calogero functional equation ()A.2n|) and the heat equation ()A.13|) . 

2.4 Isomonodromic deformations and integrable systems 

We can consider the isomonodromy preserving equations as a deformation ( Whitham quantiza- 
tion) of integrable equations [3?'. The level k plays the role of the deformation parameter. Here 
we investigate the particular example - the integrable limit of the vector generalization of PVI 

Introduce the independent time t as t = tq + Kt for k — ^ and some fixed tq . It means that 
t plays the role of a local coordinate in an neighborhood of the point tq in the moduli space A4 
of elliptic curves. It follows from (|2.17j) that the limiting curve is St-q = {z, z}. In this limit we 
come to the equations of motion of CM 1)2. llj) - 1)2. 13|) : 

dfUj = Vj , 

dtVn = - '^PjkPkjdu„E2{Uj - Un, Tq) , (2.26) 

dtp = 2[Ju(ro) • p,p] . 

The linear problem for this system is obtained from the linear problem for the Isomonodromy 
problem (|2.23|) . 1)2.24(1 by the analog of the quasi-classical limit in Quantum Mechanics. Repre- 
sent the Baker- Akhiezer function in the WKB form 

5(0) 

^' = $exp( (2.27) 

K 

Substitute (tOTl) in the hnear system (|22S1), (12:2^ . If d^S^^^ = and dtS^^'^ = 0, then the 
terms of order vanish. In the quasi-classical limit we put 95(0) ^ X. In the zero order 
approximation we come to the linear system for CM 

i. (A + L(o)(z,ro))y = 0, 
a. dzY = 0, 

iii. (at + M(o)(z,ro))y = 0. 

The consistency condition of this linear system 

dtL^^\z) - [l(o)(z),m(o)(z)] =0, 

is equivalent to the Calogero-Moser equations (|2.2H|) 
The Baker-Akhiezer function Y takes the form 

y = $6 9^0 

The same quasi-classical limit can be applied for the monodromy preserving equations that 
will be considered in next Section. 
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3 Isomonodromic deformations and Elliptic Top 

3.1 Euler- Arnold top on SL(iV,C). 

Let S e sl(iV, C)*. Expand it in the basis (lR4l) S = ^ ^^(2) ^aT^. 

According with HB.6|) the Lie-Poisson brackets on sl(A^, C)* assume the form 

{Sa, Sp}i = C(a, P)Sa+i3 ■ 

The Lie-Poisson brackets are degenerated on g* = sl(A^, C)* and their symplectic leaves are 
coadjoint orbits of SL(A^, C). To descend to a particular coadjoint orbit O one should fix the 
values of the Casimirs for the linear bracket. The phase space is a coadjoint orbit 

= {S G 0* I S = gSog-' , g G SL(iV, C) , So e Q*} . (3.1) 

The Hamiltonian of the Euler- Arnold top has a special form. It is a quadratic functional on 

0* 

H = -hr{S,JiS)), SG0% 
where J is an invertible symmetric operator 

J : 0* ^ . 

It is called the inverse inertia tensor. The equations of motion assume the form 

9tS = {J(S),S}i = [J(S),S]. (3.2) 

3.2 Non-autonomous Elliptic top (NAET). 

Consider a special form of the inverse inertia tensor. Let 



E2{a) = £^2( "^ ^"^"^ Ir) , a = {ai, 02) e 



where zj^^ is defined by ()B.3|) . Define as above the diagonal action 

J(r,S) : Sa^E2{a)Sa, J(t, S) = J(r) • S . 
Then the Hamiltonian of ET assumes the form 

^iVAET(s,^) = Zitr(SJ(r).S) = -i J2 ^7^2(7)^- 



The equation of motion of NAET is similar to (|3.2|) 

^a^S = [J(r) -8,8], ( KdrSa = S^Sa-^E2{j)C{a, 7) ) . (3.3) 



As in the previous case one can consider the limit k — > 0, (r — > tq) to the integrable elliptic 

top 

dtS = [3{ro)-S,S]. (3.4) 
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3.3 Lax representation 

3.3.1 Flat bundles of degree one 

Let be a flat bundle over the deformed elliptic curve S,- of rank N and degree 1 with the 



connections 



where L^^\w , w , t) , L^^\w,w,t) are meromorphic maps of T,r in sl(A^, C). 

For generic flat bundles of degree one the connections can be chosen in the form 

1. L^^^ = 0. From the flatness one has 

= . 

2. The Lax matrices satisfy the quasi-periodicity conditions 

L«(u; + 1) =QL«HQ-i, 
lW(u' + t) = AL^HA-i + ^ , 

T 

A{w,t) = -eN{-w - -)A 

for Q,A ()B.1|) . ()B.2|) . It means that there are no moduli parameters for E^j. 

3. -L^^) has a simple pole at t/; = and all degrees of freedom come from the residue 

Res\u,=o L'-^Xw) = S. 
The Lax matrix is flxed by these conditions: 
Lemma 3.1 The connection assumes the form 

L^^Hw) = -^d^ln^iw;T)Id+ 5„^„Hr„. (3.6) 

where ip{a, w) is defined by llB.l(jl\) . and Tq, are the basis elements iB.4[ )- 
Fixing the connections we come from H3.5|) to the linear system 

i. (Ka^ + L«H)^ = 0, . . 

As above, the independence of the monodromy of (|3.7|) means that the Baker- Akhiezer vector 
satisfies the additional linear equation 

Hi. (kS^ + M(^))^' = 0. (3.8) 



Lemma 3.2 The equation of motion of the non- autonomous top is the monodromy pre- 

serving equation for ^cl. 7| ) with the Lax representation 

drL^^^ - 9^m(i) + - [m(i) , = , (3.9) 

K 

where L^^^ is defined by \8. (^) . 

M^^^ = -^dr\n^{w-T)Id+ S^f^{w)T^, 



and f-yiw) is defined by IIB.11\) . 
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The proof of the equivalence of (|3.3|) and (|3.9|) is based on the addition formula (|A.2fl|l and the 
heat equation (|A.13|) as in the case of CM. 

In the quasi-classical limit k ^ we come to the integrable top on SL{N, C) (|S31). 

4 Symplectic Hecke correspondence 

We construct here a map from the phase space of CM (|2.7|) to the phase space of ET (|3.1j) 

:n<'^^^n^^, ((v,u,p)^S), (4.1) 
such that H"^ is the symplectic map 

H+*cu(S) = w(v,u,p). 

To construct it we define the map of the sheaves of sections r(V^°'*) — > r(V^^'*) such that it 
is an isomorphism on the complement to w and it has one-dimensional cokernel at w € S,- : 

r(i4')) ^ r(y«) ^ cu ^ o . 

It is the so-called upper modification of the bundle E^^^ at the point w. 
On the complement to the point w consider the map 

r(F«)^r(Fi°)), 

such that =Id. It defines the lower modification at the point w. 

In general case the modifications lead to the Hecke correspondence between the moduli spaces 
of the holomorphic bundles of degree k and k + 1 

j^k _^ j^k+i _ ^4 2) 

The modifications can be lifted to the Higgs bundle. They act as singular gauge transformations 
on the Lax matrices and provide the symplectomorphisms between Hitchin systems {symplectic 
Hecke correspondence). 

The modifications on the Higgs bundles can be applied for the monodromy preserving equa- 
tion as well. The action of the upper modification on the Lax matrices has the form 

L(^) = + . (4.3) 

This form of transformation provides its symplectic action. 

Consider in details the map (|4.1|) . According with its definition the upper modification H^(z) 
is characterized by the following properties: 
• Quasi-periodicity: 

H+(z + l,r) = -QxH+(z,r), (4.4) 
H+(z + r,r) = A(z,r) X H+(z,r) X diag(e(uj)), Hz,t) = -eN{-z - (4.5) 

•• ^"'"(z) has a simple pole at z = 0. Let rj = (rj^i, . . . ,rj^jv) be an eigen-vector of the matrix 
p € O (E31), pr, = pOf.. Then Res{E),=ori = 0. ' 

The former condition provides that the quasi-periods of the transformed Lax matrix corresponds 
to the bundle of degree one. The latter condition implies that L^^^ has only a simple pole at 
z = 0. The residue at the pole is identified with S. 
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First, we construct {N x N)- matrix u; r) that satisfies (|4.4j) and (|4.5|1 but has a special 
one-dimensional kernel: 



J 1 

iij[z,u;T) =U {z-Nuj,NT), 

2 J 

(z, r) is the theta function with a characteristic 



where ^ 



a 
b 



(4.6) 



a 
6 



It can be proved that the kernel of H at z = is generated by the following column- vector : 



(-1)' H ^Uk-Uj,T)\, l = l,2,---,N. 

j<k;j,k^l 



Then the matrix E{z, u, rj), (rj = (ri^j, . . . , r7v,j)) assumes the form 



H+(Z,U, Ti 



diag 



(-ly 



n ^i^k 

j<k;j,k^l 



Uj,T) 



(4.7) 



The transformation (|^ with S (|I?7j) leads to the map 7?.^*^ TZ^'^ (|^ . 

For the spinless CM, defined by the coupling constant v"^, this transformation leads to the 
degenerate orbit of the NAET. 

Note that the equation for the spin variables of CM H2.13() reminds the equation of motion 
for the NAET with the coordinate-dependent operator Ju (|2.14j) . The only difference is the 
structure of the phase spaces n^^^ ^n\, and n^'^ (inU). The upper modification H+ carries out 
the pass from 'Rp^ to TZ^^ . It depends only on the part of variables on TZ'-^^^ , namely on u and 
p through the eigenvector r^. For the rank two bundles it is possible to write down the explicit 
dependence S(u, v, u). We postpone this example to the general case of PVI in Sect. 6. 



5 Elliptic Garnier-Gaudin models 
5.1 General construction 

Here we present another way to construct the elliptic form of PVI (|1.4|) that does not appeal to 
the transformation ()6.3() of its rational form The advantage of this approach is a simple 

derivation of the linear system and the corresponding Lax operator. 

The CM and ET hierarchies are particular case of the following construction. The main 
objects are the Lax matrices and the corresponding integrable systems. Let A = {za, a = l...m} 
be the divisor of marked points on T,^ and V^(St-\A) is a holomorphic vector bundle of rank N 
and degree k over S,- with quasi-parabolic structure at the marked points [HJI . This data allows 
us to define an elliptic analog of the Garnier-Gaudin system. 

The moduli space A^C'I^I^) of holomorphic vector bundles of degree k and rank N over S,- 
is parameterized by d elements u = {ui, . . . , Ud), uj £ S,- and 

d = dim(M(°l'=l^)) = g.c.d.(A^, k) . 
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The bundles V^{T,t-\A) and V^^^(St-\A) are isomorphic. The isomorphisms is provided by 
twisting by a hnear bundle Thus, in the definition of k is defined up {jnodN). 

The transition matrices for a generic holomorphic vector bundle can be chosen as follows. For 
d = g.c.d.(A^, k) define the diagonal matrix U containing I = N/d blocks (ui, . . . , u^) 

U = diag ((til, . . . ,Md), . . . , (ui, . . . . 

Then the transition matrices corresponding to the two basic non-contractible cycles take the 
form ^ 

gi = Q , Qr = —eN{—kz — 2^t) exp(27r?U)A'^ . 

The quasi-parabolic structure means the fixing a flag structure Fla at z = Za- One can act on the 
set of flags by the group Dk = Centr{gi, gr). It is easy to see that Dk is a subgroup of diagonal 
matrices (D^ C D c Sh(N,C)) and dimZ^^ = d. The moduli space of the quasi-parabolic 
bundles is defined as 

TW^*"!^'!^) X {UaFla)/Dk. 

Then 

m m 

dimX^'"!^'^) = g.c.d.(iV, k)+Y^ dim Fla - g-c.d.(iV, k) =Y^ dim Fla ■ 

a=l a=l 

In what follows we assume that all fiag varieties are non-degenerate. Thereby 

dimA^^"^!^'!^) = ^mN{N - 1) . 

The elliptic Garnier-Gaudin models EGG{m\k\N) are particular examples of Hitchin systems 
[13j related to holomorphic vector bundles of degree k and rank N. Particular cases of the elliptic 
Garnier-Gaudin models were considered in Ref. ^Hl (k = 0) and in Ref. [201123 (k = 

The Lax matrix L('"l^'l^)(z)dz is a meromorphic section of the bundle Endy^®r2*-^''^^(ST-\A) 
with simple poles at {za} G A and fixed residues 

i?esLH^1^)(2j = S'^ G s\*{N,C) ■ (5.1) 

More exactly, we assume that S" G O"" is an element of a non-degenerate coadjoint orbit. 
The Lax matrix L^™!*'!^) (21)^2 plays the role of the quasi-parabolic Higgs field in the Hitchin 
construction ^. 

The Higgs field respects the lattice action Z rZ: 

LMk\N)^z + 1) = g-^L^--\m)^z)g^ , Lim\k\N)(^^ ^ ^) ^ g^^ L^"'\>'\^\z)gr . (5.2) 

It follows from 1)5. 2|) that the diagonal part of L^"^^''\^\z) is a double-periodic meromor- 
phic function. It depends on the additional variables {vi, . . . ,Vk)- The conditions (|5.1j) . (|5.2|1 
determine the Lax matrix L^"^^^\^\z). Below we consider particular appropriated cases. 

The whole phase space 

j^Mk\N) ^r,.*^(o\k\N) ^ = {v,u;p\...,p™} 

is equipped with the Poisson brackets 

{vj,Uk} = Sjk , {Pu , P% } = ^'"^ {Ptjij - Pji^ik ) • 

^In what follows we omit the differential dz in the notation of the Lax matrices and work with L{z) instead of 
L{z)dz 
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In the case when dim(X(0|'=l^)) / one can go further and consider the symplectic quotient 
with respect to the defined above group 

^H^^IA^) ^^Hfc|A^)//l),. (5.3) 
The dimension of the reduced space 

m 

dim7^^7J'''^^ = ^ dimO'" = mN{N - 1) 

a=l 

is independent on the degree of the bundles. In fact, ah the spaces '^j,^^'^'' for 



k = 0,...,A^ — 1 are symplectomorphic [22] . The symplectomorphisms is provided by the 
symplectic Hecke correspondence placed between the Higgs bundles of degree k and k + 1 ()4.2() . 

The commuting Hamiltonians of EGG{m\k\N) can be extracted from the expansion over 
the basis of the double-periodic functions Ei{z — Zg) (|A.2|) . HA.4|) . (|A.5|) 

\i{L^'^\^\^\z)y - J' Hf^Ei{z - Za) , {j = 1) 



j I K + T.a ELi Hf^E.,{z -Za), ij = 2,...,N), 



Oj Z—ia Z-/i=l ij 

where H'^j = 0. Equivalently, the commuting Hamiltonians can be defined by the integral 
representations ^Sl 

Hfj = jJ^ tr(LH'=l^)(^)y^^.. 

Here = lJ.ij{z, z) are (1 — j, 1) differentials on St-. They are chosen to be dual to the basis of 
the elliptic functions 

- / f^tjEiiz - Zb) = 6f!' . 
3 J S7- 

They have the following form j2Ul I21j . For a marked point Za define a smooth function Xa{z, z) 
as for z = (t2J6l) Then for a 7^ 0, i 7^ 

^i% = ^^{z-Za^^^Xa{z,z) 



and 



° - Kd{-Z-z){l-Y,Xa{z,-z)). 



T -T 



The higher order differentials ;Ug -, j > 2 are related to the so-called W-structures on S,- |2^ and 



will not be discussed here. 

The commuting Hamiltonians generate the hierarchy of the Lax equations 

5(„,,)L(-I^I^)(.)) = [L^-\^\^\z)l m(-I^P(.)] , = 9,^^,,,^, . (5.4) 

Here the matrices M^'^^^j^^'^ have the same quasi-periodicity properties as L^™'!'^!^) and satisfy 
the equation 

The Lax equations (|5.4jl as well as the form of the M-matrices follow from the equations of 
motion for the Hitchin systems. Their derivation can be found, for example, in j28j . 
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5.2 Involution of the Higgs bundles 

Let ? be an involution {q^ = Id) acting on the space of sections of the Higgs bundle. Then we 
have decomposition 

r(Endv^) = r+(Endy^) r-(Endy^) . 

Let 

Assume that the involution preserves the hierarchy 

and we choose the invariant subset {M^^^^^''^}. Then the constraint 

= (5.5) 

is consistent with the involutions (|5.4j) 

Consider the case of the special divisor 

A = (wo = , a;„ , a = 1, 2, 3) . (5.6) 

Then the transformation z — > —z preserves S,- \ A. It can be accompanied with an involution 
of sl(A^, C) to generate <j. In next subsections we consider two examples of the rank two bundles 
with the involutions of such type. 

5.3 Degree zero bundles 

We have derived the space A^C^I'^I^) in 2.3.2. Remind that it is described by N parameters 
u = (ui, . . . , Mat), G St-, such that ui + ... + mat = and the multiplicators are: 

gi=Idiy, (^T- = e(— u) = diag(e(— Ml), e(— Mat)) , (e(M) = exp 27rzn) . (5.7) 

These conditions fix the Lax matrix of EGG{m\0\N) up to a diagonal matrix diag v, 
V = (f 1, . . . , uat). It follows from (|5.2|) and (|5.7|) that the Lax matrix assumes the form 

m m 

L\-^m\z) = 6,,{V, + Y^PuM^ - ^a)) + (1 - 5i,)Y,P%<t^{^^ - ^ - ^a) ■ (5.8) 
a=l a=l 

The particular important case corresponding to a single marked point was considered in Section 
2. 

Now consider the special form of (|5.8() with A ()5.6|) and N = 2. The Lax matrix 1)5. 8(1 can 
be gauge transform to the form 

3 

LW0|2)(^) ^ (^ + ^^a^^(^ + ^jy3+ (5.9) 
a=0 
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a=0 

Define the involution as 

?(0)(M(4|0|2)(^)) = c7iM(4|0|2)(_^)ai . 



(5.10) 



Proposition 5.1 ( j43| ) The invariant Lax matrix L'^^'^^'^ ^ with respect to <j^'^^ assumes the form 

l(^|0|2) + (z) = L^^(z), 

3 



L 



CI 



V 
-V 



a=0 



z:;a99a(2u,Z + Wa) 

DaiPa{-2u, Z + Ua) 



, (5.11) 



w/iere = ^trCp")^; V'aC^^^) + i^a) = e(-2n5rWa)0(2u, ^ + Wa). 
Proof. 

The projection of the Lax matrix on its anti-invariant part should vanish. It follows from (|A.9|) 
and (|A.17I) that 

e{2udrUJa)4>{2u, Z + OL>a) — > - e(-2u5ru;a)'/'(-2ti, Z + Wa) • 
Since aia^ai = —a^, a\a^a\ = (T_ the vanishing of L^^l*^'^) ~(z) implies 

P3 = 0, P^=P^, (a = 0,...,3). (5.12) 
Thus one can put p\= p°L = Va and the invariant part L*^^''^'^) ^(z) coincides with (|5.11|) . □ 

Remark 5.1 Remind that upon the involution the phase space i5.S\) is the symplectic quotient 
'^red'^'* ~ 7^*-^''''^V/^2- The Conditions 115. 1^) imply the moment constraint X^a^'a ~ ^ '^"'^ 
gauge fixing X^^p^ = ^aP- ^'^'^ thereby provide the pass to the symplectic quotient. 



According with Ref. j43j the new Hamiltonian can be read off from the decomposition 



a=0 



(5.13) 



where 



and 



a=0 







/ 1 


1 


1 


1 \ 








1 


1 


1 


- 1 


- 1 








~ 2 


1 


- 1 


1 


- 1 






\ ^3 I 




V 1 


- 1 


- 1 


1 ) 




V ^3 / 



(5.14) 



The canonical bracket {v^u} = 1 and the Hamiltonian H'~"^ yields the BCi CI system (|1.5j) . 
Proposition 5.2 ( |43| ) There exists the matrix M^l*^l^+ = M'-^^ , 





i>a^a{-'^U, Z + UJa) 



a=0 

SMc/i ^/iat t/ie Lax equation is equivalent to the BCi CI equation \1.5\) 



^The difference of the ? action on L and M is due to L G n(^'"'(E^) while M € r2'°^(Er) 
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5.4 Degree one bundles 

The Lax matrices in this case was considered in pHI- The multiphcators of a bundle are: 
91 = -Q, 5r = -Cat (^-^ - 2;^ A, eAr = exp^ 

where Q and A are the matrices ()B.1|) . (|B.2|) . There are no moduh parameters in this case. The 
phase space is a direct sum of the coadjoint orbits x ... x and 7^(''"li|^) = TZ^^'^^^\ In 
the basis {T^}, (a G ^^^) <|B.4|1 the Poisson brackets are 

{5S,4}i = C(a,/3)<5'^^5S+/3- (5.15) 
The Lax matrix is completely fixed by the quasi-periodicity conditions 

m— 1 

^Hi|A^)(^) = ^ ^ r,5>,(z-z,), ^^{z) = e{zdruMz,u^), uj^ = ^l±pl, (5.16) 
Consider the case N = 2 and A 1)5.61) 

3 3 

^(4|l|2)(^) = ^ SM^ - ^a)c7a • (5.17) 
a=0 0=1 

Define the involution 



«(L4|i|2(z)) = -L4|i|2(_z), ?«(M4|i|2(z)) =m4|i|2(-z). (5.18) 



Proposition 5.3 The invariant Lax matrix L^^'^'^^"^ with A J ,'5. 6|) assumes the form 

lWi|2) + (z) = L^^«(z). 

a=l a=l ^ I \ / / 

where Sa = S^, i'a = S"^ = const and 



Proof. 

It follows from ifOT^ and (fRTsl) that 



L^'\'\'^~ = i{L^'\'\'\z) - l^iii2(-.)) = y: SM 



(5.19) 



(5.20) 



If 2^{4|i|2) _ then S"^ = for a ^ a. Fixing the Casimir functions before the involution as 
^^(5^)2 = v"^ we conclude that 5^ = On the other hand 

3 

^(4|i|2)+ ^ ^ (5°(^„(z) + 5>„(z - a;„)) (7„ . 

Q = l 
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and we come to the statement. □ 

Again one can construct Hamiltonians from the Lax matrix as in H5.13() . Then the 

Hamiltonian 

i^^^^ = -Ytr(S,(J(r).S + z.')) (5.21) 

with the H5.15() Lie-Poisson brackets leads to the ZVG equation (|1.9() . 

It is worthwhile to mention that we deal with three types of constants v and v' . The 
latter two are expressed through the former by H5.14() and H5.20() . 

Proposition 5.4 There exists the matrix M^^l^l^^"^ = M^^'^, 

M^^^iz) = J2 [-Sa^a{z){Ei{z + uja) " E^iuj^)) + v^z - uj^)Ei{z)\ a« = (5.22) 

a 

= E -S.^liiM^a. + E,iz)L^y^iz) . 

such that the Lax equation is equivalent to the autonomous SL(2, C) ZVG \6.1)) . 
The proof of this statement essentially is the same as Proposition 5.2. 



6 Non- autonomous systems 

Consider the non-autonomous systems NAZVG H1.2|) . and EPVI 1)1.4(1 

drS = S^J{T,S) + S^u', (6.1) 

-l^ = -Y.ylp'{u + u:a,T). (6.2) 

a=0 

We remind the relation between the last equation and PVI ((1.11) . established in 23.^. Assume 
that the parameters of the equations are related as follows 

fo = a , vi = -13, V2 = l , 1^3 = ^ - • 
Then the substitution in (|1.1|) 

(n, r) ^ X = , t = \ ea = E2{uJa) (6.3) 

\ 62-61 62 - ei / 



brings it in the form (|6.2|) . 

The Lax matrices for the non-autonomous systems can be obtained from the Lax matrices 
for integrable systems, constructed in previous section, by replacing the coordinates z ^ w 
(|2.17() . In this way we come to the following two statements. 

Proposition 6.1 ( |43| ) The equation of motion of EPVI h6. has the Lax form 

drL^^ - d^M^^ + -[M^^,L^^] = (6.4) 

with 

L^^ = P + X , P = diag(z;, -v) , (6.5) 
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Xi2{u,w) = '^DaiPa{2u,W + UJa) , X2l(n, w) = Xi2{-U,w) , 
a 

M^' = Yjk, ij^k,j,k = l,2) (6.6) 

Yi2iu,w) = '^Daip'a{2u,W + iOa) , Y2l{u,w) = Yi2i-U,w) . 
a 

Proposition 6.2 The equation of motion of NAZVG has the Lax form with 

j^NAZVG ^ _^Q^ ^ + vo^^^iw - u^o)W . (6.7) 

a 

^NAZVG ^ _^Q^ + y _^ Ji{w)M^)y.sH ^^ ^ E,{w)L^^^^^{n = 0) . (6.8) 

The proofs of the both statements repeat the autonomous cases. The only new ingredient is the 
heat equation ()A.13|) . that should be used. 

Now we establish the interrelation between these systems. It is provided by the same modi- 
fication transform as above. Namely 

lNAZvg^^^ u;) = H+Aca(H+)-i + ^ 

where H"*" is defined by (|4.7j) . 

The upper modification leads to the following relations between the coordinates of the phase 
spaces 

Q - „. g2(0) e2(2n) K 02(0)^2(2") I 
'^1 ~ '^i?'(0) i?{2«) 2 ,?'(0) i?(2«) 

~ Bl{0) ez(2u)ei{2u) ~ el{2u) ~ 02(0) 6>2(2m)6>4(2m) ~ 62(0) 6>2(2m)6>3(2m) 
^093 (0)6*4(0) W^2^) ^ '^1 i?^(2u) '■^20^ Jj7(2^^) \- I^Sg^ W^(2u) ' 

„-o _ „,e3(o) e3(2«) ^e3(o)eg2«)_ 

'^2 — t'^^T^ ^(2«) "T 2 iFtSy i9(2m) /p qn 

eg(o) 62 (2^)04 (2m) _ 63(0) e3(2u)e2(2u) _ ~ 03(0) 63 (2^)04 (2m) _ ~ ei(2M) 

6*2 (0)6*4 (0) 5?^(2^I) '^1 6*2 (0) ;?^(2^I) (q) 5j^(2;I) '^3 ^2(2^) , 
Q - -,. ^4(0) 6*4(2m) k 6^4(0) e'^{2u) 
'-'3 ~ "^iFtoy i?(2m) 2 iFtSy i?(2m) 

- e|(0) e2(2M)e3(2n) - 64(0) e2(2n)e4(2M) , e|(2n) ~ 64(0) 04 (2^)03 (2«) 

^0 6*2(0)6*3(0) 5?^(2^I) ^ '■^16*2(0) ^&^(2u) ^ '^'^W{Iu} ^35)3(0) ^^^(III) ■ 

The proof is based on the direct usage of ()B.19|) . ()B.2n|) and Riemann identities for theta 
functions ^5] . 

7 Quadratic brackets and NAET 
7.1 r-matrix structure 

The classical r-matrix is the quasi-periodic map from to End(£:^^VEnd(^JJ^) g] 

r{w) = ^ip.,{w)T^®T..,, (7.10) 
7 

where 93^ is defined by HB.10|) . It satisfies the classical Yang-Baxter equation 

[^(12) _ ^')^ [w)] + [r(^2) _ ^(23) (^/)] 
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+ [r^^^\w),r^'^^\w')]=0. 

By means of the r-matrix one can define the Unear brackets. Let L^^\z) be the Lax matrix (|3.6|1 
for K = 

7 

Proposition 7.1 ( jl7| ) The Lie-Poisson brackets onsl{N,C) 

{Sa,Sf3}i = C{a, l3)Sa+/3 

in terms of the Lax operator L^^\S, z) are equivalent to the following relation for the Lax operator 

{L^l^ {w), {w')}i = [r{w - w'), L(i) {w) ®Id + Id® {w')] 

L^^^=L^'^®Id, L^^^=Id®L^^\ 
The proof is based on the Fay three-section formula (|A.19|) . 



7.2 Quadratic Poisson algebra 



In addition to the A^^ — 1 variables S = a G ^j^^} introduce a new variable Sq and the 



GL(A^, C)-valued Lax operator 

L = -SoId + L^^\S,w). 

It satisfies the classical exchange algebra: 

{Li{w),L2{w')}2 = [r{w-w'),Li{w)(^L2{w')]. (7.11) 

These brackets are Poisson, since the Jacobi identity is provided by the classical Yang-Baxter 
equation. 



Proposition 7.2 The quadratic Poisson algebra i7.11[ ) in the coordinates {Sq, S) takes the form 
{5„, Soh = E S^^^S^iE^C-^) - E2i^^^^^^l±p^))Cia,,) , (7.12) 

{Sa,Sf3}2 = SoSa+f3C{a,(3) + ^ Sa-ySi3+^f{a,(3,j)C{j,a-(3), 

where 

f(a,A7) = E,0^) +E,i ^^+^^-^^ +lf' + - "^)- )- 

f Pi+li + (/32 + 72)t- ^ , „ ai -71 + (a2 -72)^ ^ 
Ei{ ^ ) + Ei{ ) . 

It is the classical Sklyanin-Feigin-Odesski (SFO) algebra jHSllHI- These brackets are extracted 
from (TtTTI) by means of (TOTl) . (fX22|l . 

Two Poisson structures are called compatible (or, form a Poisson pair) if their linear com- 
binations are Poisson structures as well. It turns out that the linear and quadratic Poisson 
brackets are compatible, namely, there exists the one-parametric family of the Poisson brackets 

{S,S}a = {S,S}2 + A{S,S}i. 
22 



In the case of the ET integrable hierarchy these compatible brackets provide the hierarchy with 
the bihamiltonian structure |^. The hierarchies of the monodromy preserving equations are 
more intricate and we do not consider here the hierarchy of NAET. However, we have the 
fohowing manifestation of the bihamiltonian structure. 

Proposition 7.3 In terms of the quadratic brackets the equation of motion of NAET 
assumes the form 

K-drSa = {So, Sa}2 ■ 

The proof follows immediately from (|7.12j) . We replace the linear brackets on quadratic and 
simultaneously the quadratic Hamiltonian H^^^"^ on the linear Sq- 

8 Reflection equation and generalized Sklyanin algebra 



In this Section we present another Hamiltonian form of NAZVG (|1.2|) . It is based on the 
quadratic Poisson brackets. The quantization of these brackets is described by quantum reflec- 
tion equation. 

8.1 Quantum reflection equation 

Let R~ be the quantum vertex R-matrix, that arises in the XYZ model. We introduce also the 
matrix i?"*" 

3 ft 



R^{z,w) =''^^ipa{z ± w)aa ® (Ta (8-1) 



a=0 



where ipi is defined by ()(].H) . Define the quantum Lax operator 

L{z) = So(l)\z)ao + Y^iSa^^iz) + D^ip^{z - c^„))a„ . (8.2) 

a 

Proposition 8.1 The Lax operator satisfies the quantum reflection equation 

R-{z,w)Li{z)R^{z,w)L2{w) = L2iw)R^ {z,w)Li{z)R- {z,w) . (8.3) 
iff its components Sa generate the associative algebra with relations 

[9a,up]=0, [P„,5a]=0, (8.4) 

i[So,Sa]+ = [Si3,S^/], (8.5) 

[Sj,So] = i^^--^[Sa,Si3]+ - 2i^{i>aPaSf3 - vpPpSo,) , (8.6) 

where 

Ka = Ei{h + UJa) — Ei{h) — EiiuJa) , Pa = — e^p{ — 2TTlUJadri^a)4'{^a + fi, —^a) ■ 
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For the proof see Appendix D. 

If all pQ, = the algebra ()8.4() - (|8.6|) coincides with the Sklyanin algebra. Therefore, it is a 
three parametric deformation of the Sklyanin algebra. 

Two elements 



a 1 



a 

belong to the center of the generalized Sklyanin algebra (|8.4() . (|8.5I) . They are the coefficients of 
the expansion of the quantum determinant 

det(L) = trP(L(z, h) ® L+{z, -h)) , 
h 

where 

a 

a 

over the basis of the elliptic functions. However, we do not know how to derive the expression 
for the quantum determinant directly from the reflection equation. 

8.2 Classical reflection equations 

Consider (|8.3j) in the limit ?i — > 0. The classical r ^-matrices are defined from the expansion 

R^{z,w) = (-j ao0ao + r^{z,w)+O{h), 

r^{z, w) = ^ tfa{z ± w)aa ® CTq, . 
a 

For the Lax operator one has 

L{z) = h'^Soao + '^{Sa(Pa{z) + i>a^a{z - Wq,))^^ + 0{h) . 
a 

Define the classical variables 

Sa ^ Sa , So ^ hSo , 

and the corresponding classical Lax operator 

3 

L{z) = Soao + '^{Sa(Pa{z) + Ua-^aiz - UJa))aa ■ 
a=l 

Then, taking into account that [Li,L2] = h{Li,L2} + 0{h'^) one finds the classical reflection 
equation in the first order of 

{Liiz),L2{w)}2 = ^[Li{z)L2{w),r'{z,w)] + (8.7) 
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-L2{w)r^{z,w)Li{z) - -Li{z)r^{z,w)L2{w) . 
by passing from the group- valued element L. For the Lie-algebraic element L 

3 

a=l 

we come to the linear brackets 

{Li{z),L2{w)}i = -^[r-{z,w),Li{z) + L2{w)] + ^[r+ {z,w), Li{z) - L2H] . (8.8) 

Proposition 8.2 The classical reflection equations \8. 7\ ) leads to the quadratic Poisson struc- 
ture on generalizing the classical Sklyanin algebra 

{Sa, Si3}2 = 2ieaf3-ySoS^ , (8.9) 

{5*0, Sa}2 = ieaf3'YSf3S^{E2{LOf3) - -^2(0;^)) + 2i£ ^fS^ S pv'^ , (8.10) 

where v' is defined as l^5.2(J\} . 
Proof. 

The only dissimilarity from the classical Sklyanin algebra is the linear term in (|8.1fl|l . It comes 
from the last term in the right hand side of (|8.6|) in the limit h ^ 0. ID 
There are two Casimir elements of the Poisson algebra (|8.9j) . (|8.10|) 

a 

C2 = 8^ + J2(pMSI + 'i^'aSa) ■ 

a 

They are the coefficients of the expansion of det(L(z) over the basis of elliptic functions. 

It is easy to see that the brackets ()8.9|) are compatible with the linear SI2 Lie-Poisson brackets 

{Sa, Sff\i = 2ie„^^5-y . 
The equation of motion of NAZVG H1.2|) is written in terms of the linear brackets: 

drSa = {H^^'^, Sa}l ■ 

The straightforward calculations shows that (|1.2|) can be written in the form 

drSa = {Hq, Sa}2 , Hq = Sq . (8.11) 

In this way for the generic form of PVI we have the same analog of the bihamiltonian property 
as in the degenerate case (Proposition 5.4). 
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8.3 Spin chain with boundaries 

Quantum reflection equation allows us to define XYZ model on a finite lattice with boundary 
conditions |31]. The Lax operator (|8.2|) can be considered as a new solution of the reflection 
equation. 

Consider a pair of matrices K^{z) with Poisson brackets 
{Kt,Ki} = [Ktiz)Ktiw),riz-w)]+Ki{w)riz + w)Kt{z)-Kt{z)riz + w)Ki{w) (8.12) 
and L^{z), i = 1...N with brackets 

{L\{z),Li{w)} = 5^^[r{z - wlL\{z)L\{w)] . (8.13) 

Then, according with [M] 

h{z) = tr \k+{z)L^{z)...L\z)K-{z) {L\-z))~^ ... (L^(-z))"^1 (8.14) 
is the generating function of commuting Hamiltonians 

{h{z),h{w)} =0. (8.15) 

Choosing L^'^^± (2) = L^'^G± £qj, j{^{z) we construct a spin chain with boundaries. But 

in (|8.7|) we have a factor ^ which comes from (|8.1|) on the quantum level. Thus, we should put 
the brackets on boundaries for L^^*^±(z) to be two times more than in ()8.9() . In other words 
i?-matrices for L^^*^±(z) should depend on the same Planck constant h as for all other U{z). 

Proposition 8.3 Spin chain involving N internal vertices L^{z) with boundaries L^^'-'^[z) is 
integrable if 

{S^,St} = 2teo,p^iSp^,iE2iu;p)-E2iu;^)) + Sfiy'^) > 

and for i,j = 1...N 

{Si Si} = 5'heo,p^SlS\{E2{u:p) - E^iuo^)) . ^ ' ' 

The nearest-neighbor interaction is described by the Hamiltonian: 

H = ln (s^S', + E S~Si{C - E2M) + {i^XSi] + 
In i^S^S+ + E S^S+{C - E2M) + ii^')+S^j + (8.18) 

N-l / \ 

E In 5^5^+^ + E S'i,S:,+HC - E2M) , 

1=1 V a / 

where C is a constant, equal to the fraction of the values of the Casimir functions for Sklyanin 
bracket 118.1', 



The proof of 1)8. 18() is similar to those one given in [7] (Chapter III, §5). Consider a special 
point E S, such that E2{zq) = = Cj and assume that this relation 

Eqi qi 

a = l 

is independent on a point of the lattice Cj = Cj yi,j = 1...N. In this case all L*(zo) are 
degenerated and have the form L^{zq) = a*(zo) x P^izo) for some vectors a*(2;o) and covectors 
P^{zq). To finish the proof we should notice that all detL*(2;) are Casimirs for 1)8. 13() and in our 
case U{z)U{-z) = 1 • detL^(z) Vi. 

Some solutions of the reflection equation corresponding to the spin chains with dynamical 
boundary conditions are constructed in Ref. |19| I15j. 
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9 Appendix 

9.1 Appendix A. Elliptic functions. 

We assume that q = exp27rfr, where r is the modular parameter of the elliptic curve E-,-. 
The basic element is the thcta function: 

d{z\T) = q^^{-lY-e{^n{n + l)T + nz) = (e = exp27rz) (A.l) 
nez 

oo 
n=l 



The Eisenstein functions 



Ei{z\t) = d,log^z\T), Ei{z\t) r~.^-2r]iz, (A.2) 



where 

oo oo 



_3_ ^ V" 1 _ 24 V(r) 



m=— oo n=— oo 

where 



■>^^=— oo n=— oo ^ ' < v / 



n>0 

is the Dedekind function. 

:j2 



1 



E2{z\t) = -d,E^{z\T) = dt\og'&{z\T), E2{z\t) ~ ^ + 2r?i . (A.4) 



z 



The highest Eisenstein functions 



(A.5) 



Relation to the Weierstrass functions 

C{z,T)=Ei{z,T) + 2r]i{T)z, (A.6) 

p(z,T)=£;2(z,T)-2r/i(r). (A.7) 
The next important function is 

(f){u, z) = (j){z,u) , (f){—u,—z) = —(f){u,z). (A. 9) 

It has a pole at z = and 

4>{u,z) = ^ + E,{u) + ^{Ef{u) - p{u)) + .... (A.IO) 

Let 

f{u,z) = du(l>{u,z). (A.ll) 
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Then 

Heat equation 

Quasi-periodicity 



f{u,z) = (j){u,z){Ei{u + z)- Ei{u)) 



dr(p{u,w) - -r—.dudw(t){u,w) = 0. 



-d{z + 1) = -^{z) , + r) = -q-2e-^^'^^{z) 

Ei{z + 1) = Ei[z) , Ei{z + r) = Ei(z) - 2m , 

^2(^ + 1) =^2(^), ^2(2 + r) =E2(^), 



(t){u, z + l) = (t>{u, z) , z + r) = e"^™0(n, . 

/(n, z + 1) = /(^x, z) , f{u, z + t) = e-2™/(u, z) - 2t:i<P{u, z) . 
The Fay three-section formula: 



(A.12) 
(A.13) 

(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 



(i){ui,Zi)(i){u2, Z2) - 4>{UI + M2, Zl)(t){u2, Z2 - Zl) - (j){ui + U2, Z2)4>{ui, Zi - Z2) = . 

Particular cases of this formula is the Calogero functional equation 

z)dy(f){v, z) - (p{v, z)du(t>{u, z) = {E2{v) - E2{u))(p{u + v,z) , 
Another important relation is 
(/>(?;, z — w)(t){ui — V, z)4>{u2 + f , u)) — — U2 — v^z — w)(j){u2 + V, z)(p{ui — v,w) ■ 

(j){ui,z)(l){u2,w)f{ui,U2,v) , 

where 

f{ui,U2,v) = C{v) - C{ui -U2-V) + C{UI -V)- C,{U2 + v) . 

One can rewrite the last function as 

'd' {<d)^{ui)^{u2)^{u2 -ui + 2v) 



f{ui,U2,v) 



■d{ui — v){}{u2 + v)'d{u2 — ui + v)'d{v) 



Theta functions with characteristics: 
For a, 5 G Q by definition: 



a 
b 



In particular, the function •& ()A.1|1 is a theta function with characteristics: 



'd{x, t) 



1/2 
1/2 



(x,r) 



Properties: 



a 
b 



{z + l,r) = e(a) 



a 
b 



(A.19) 
(A.20) 
(A.21) 

(A.22) 
(A.23) 



(A.24) 
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a 
b 



{z + a'r, r) = e ( -a'^^ - a'{z + b) ' 





a + a' 




b 



a + j 
b 



The following notations are used: 9 



a/2 
6/2 



and ■!? = ^11. 



9.2 Appendix B. Lie algebra sl(A^, C) and elliptic functions 

Introduce the notation 



and two matrices 



2TTi 

eN{z) = exp(— z) 



Q = diag(eAr(l), . . . , eAr(m), . . . , 1) 
/ 1 ••• \ 



A 



1 




\ 1 







1 
/ 



(B.l) 



(B.2) 



Let 



^S^^ = {Z/NZ e Z/NZ) , zj^^ = zg^ \ (0, 0) 



(B.3) 



be the two-dimensional lattice of order A^^ and N'^ — 1 correspondingly. The matrices Q'^^A"''^ 

E Z 

z(2) 



(2) ~ (2) 

a = (01,02) € Z]y generate a basis in the group GL(A^, C), while Q'^^A°'^, a = (01,02) E 



generate a basis in the Lie algebra sl(A^, C). Consider the projective representation of Zjy in 
GL(iV,C) 

a^T, = ^e^(^)Q-^A»% (B.4) 



Tan 



N ax b.^ , , , , , 

TT-^NK 7;—)Ta+b , (a X 6 = 0162 - 0261; 

2tti 2 

2(^Ci) 



(B.5) 



Here y1^n{—^^) is a non-trivial two-cocycle in H'^{Zj^ ,7j2n)- It follows from ()B.5|) that 



where C(a,/5) = ■^sin-^(a x (3) are the structure constants of sl(A^, C). 
Introduce the following constants on Z^^^: 



(B.6) 



^9(7) = i}{ 



71 + 72T ^ 



^i(7) = £;i( ), E2{j) = E2{ 



N 



N 



and the quasi-periodic functions on 



(j)^{z) = (f>{ 



71 + 72T 



z) 



V9^(z) = eAr(72z)0^(2;) 



(B.7) 
(B.8) 

(B.9) 
(B.IO) 
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a 


(1,0) 


(0,1) 


(1,1) 






0-1 


0-2 


half-periods 


a;i = i 


^2 = 5 


— — 


(/7q,(z) 


62(0)61(2) 


64(2)6(^(0) 
64(0)61(2) 


63(^)6i(0) 
63(0)61(2) 



(B.ll) 

m = er,{j2z)M^KE,{:^l±^ + z)- E.O^^^)) . (B.12) 
It follows from ()A.8|) that 

tp^{z + 1) = eN{l2)'f'y{z) , ip^{z + t) = eN{-li)^-i{z) ■ (B.13) 

f^{z + l) = eN{l2)f^{z), f.,{z + t) = eN{--yi)f-y{z) -2m(p^{z) . (B.14) 

SL(2,C) case 

For SL(2, C) instead of Tq we use the basis of sigma-matrices 

ao = Id, ai = inTo^i , (J2 = ivrTi^i , 0-3 = -iirTifi , (B.15) 

{o"a} = {o"o, o-q} , (a = 0, a) , (a = 1, 2, 3) 
cji - i(J2 ai + ia2 

+ 2 ' 2 

The standard theta-functions with the characteristics are 



"0,0 — fa , Pl,o — (72 , t/Q,! — P4 , — Oi 



2^x, 2T)e,iy, 2r) = ^(2+^, r)02(^, r) + e2(^, r)^(^, r) 
203(x, 2r)02(y, 2r) = ^(2^, r)^(2^, r) + ^2(2+^, r)^2(^, r) 

2es{x, 2T)es{y, 2t) = e^{^,T)e^{^,T) + e^{^,T)e^{^-^,T) 

2e2{x, 2T)e2{y, 2t) = 03(2+^, r)03(^, r) - ^^(2+^, r)04(^, r) 



(B.16) 



/ T?'(0) \^ 

V5Q(^;)'/'a(^; - Wq) = -ei(-u;Q,a,-u;Q) -- — - (B.17) 

iPa{-Z - Uh) = -ifaiz - u;a)Sab + (1 " Sab)'^aiz - UJb) . (B.18) 

Formulae with doubled modular parameter: 

O^x, T)e3{y, r) + e^iy, T)e^{x, t) = 2ei{x + y, 2r)^4(x - y, 2t) 
^4(x, r)03(y, r) - ^4(y, t)03(x, r) = 2i?(x + y, 2r)i9(x - y, 2t) 
^3(x,r)e3(y,T) + ^4(y,r)04(x,r) =2^3(2: + y,2r)^3(a:-y, 2r) ^ ' ^ 

^3(3^, r)e^{y, r) - e^iy, T)9i{x, t) = 2e2{x + y, 2r)^2(a: - y, 2r) 



(B.20) 
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9.3 Appendix C. Deformed elliptic functions 

ip1{z) = eAr(a2z)0( "^ ^"^"^ +r],z), a G r] eT^r- 

It follows from ()B.13|) that 9?a(-z) is well defined on Z^"*: 

fa+ci^) = 'Pai^) > for ci,2 G Z mod iV . 
+ 1) = eN{a2)ip'l{z) , ^l{z + r) = e^(-ai - Nri)^l{z) . 

The following formulae can be proved directly by checking the structure of poles 
periodic properties: 

'/'('"^>^)V'a(^^ - ui) + 0(-'u;,?7)v9^(z + = (j}{z,r]){ipaiz - w) + ifaiz + w)) 

(f)^{z - w)(l)^{z + w) + ip'^{z - w)(pa{z + W)+ ip^^{z - w)ipp{z + w) + 

iplf{z - w)lp'\{z + w) = 2(/)2'?(z - w)(l)'"{2r]) + 2(^2»)(^ _^ 'u;)0-«'(2r/) 

<j)^{z - w)(P{z + w) + Lplt{z - w)Lplt{z + w) - i/l^{z - w)Lp^p{z + w)- 
ip\{z - w)ip]{z + w) = IcjP^'^iz - u;X(2r/) + 2,/.2^(z + u')(/j-"'(2j?) 

(pp{z - w)(j)^{z + w)- (f)^{z - w)ipp{z + w)- ipk{z - w)(fy{z + w)+ 
^"^{z - w)ipl{z + w) = 2^f{z - wM{27]) - 2^f{z + w)^-^{2l^) 

(/'^{z - w)ip\{z + w) + iplf{z - w)i/'p{z + w)- (f^{z - w)4>'^{z + w)- 
(P^'iz - w)ipl{z + w) = 2ipl^{z - w)ip^{2r]) + 2ipl''{z + w)ip~''{2rj) 

ifliz + w)^^zW{w) + ^l{z + w)<j)''{z)^']{w) = 

^"(2 + W)^1{Z)^1{W) + (^^(Z + W)^1{Z)^1{W) 

^liz - w)ip:^iz)cP^{w) -ifliz- w)cP^{z)ip:^iw) = 
-4>^{z - w)(pl{z)ip^^{w) + - w)'/I^{z)lpI{w) 

{Ei{7] + p) + Ei{rj - p) - Ei{7] + a)-Ei{ri- a))x 
((^^(z + w)ip'!I,iz)(j)'^{w) - cP'^iz + u')0^(z)(^IJ(u;)) = 
(^i(r/ + 7) + Ei(r/ - 7) - 2^i(r/))x 

{-(faiz + + V9^(Z + w)ip'^p{z)ip'l{w)) 

(j}'^{w){-Lpl{z - UJc^l^^piw - LOp)(f%{z -W)+ l/l^{z - UJi3)ip'l{w - LOa)'/ii{z - w)) = 
-(p~'^{w){ipl{z - UJa)(p'p{w - LOp)ipl{z + W) - Lp^p{z - U;p)ip2,{w - Ua)'/p{z + w)) 

(p^{w + (jJa){ip^{z — UJ^)<t)'^{w)Lpp{z — w) — <j)'^{z)iplf{w — Ld^)ip^{z — w)) = 
(j)^{-W + UJa){'P^{z - Ul^)<f>'^{w)(p'^p{z + W) + <f>'^{z)ip^{w - LO^)Lpl(z + w)) 
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9.4 Appendix D. Comments to Proof of Reflection Equation 

Here we give some comments on the proof of the Proposition 3.2. 

A direct substitution of (|8.1ll8.2j) into (|8.3|) yields three types expressions proportional to 
a 1, 1 ® fj and a ® a. Consider, for example, 1 ^ which contains the additional constants. 
By the usage of ()('.5II(].8|) it simplifies to 

[S^, So] {2^'^{h)i^^iz)(j)\w)ip^{z -w)- ct)\z)^^{w)c^\z - w)) + 
2cl)-'"{h){ip^(z)(t)^{w)^^{z + w)- (t)\z)ip^{w)ct>\z + w))) + 
^o] {2<t,'"{h){^^{z - uj^)(t>''{w)^^{z -w)- <i)^{z)^^{w - uj^)(t)\z - w)) + 
20-"'(/i)((^^(z - uj.i)(P\w)^'}^{z + w)- (l)\z)ip^{w - uj^)(l)\z + w))) + 

AiVaSp (^^%{z - LOa)^^{w){Lp^{z - w)(t){h, w) +ip^{z + w)(l){h, -w))- 

Lp^p{z)Lp%{w - - w)(l){h, w) + 93^(2 + w)<t){h, -w))^ + 

'iii'pSa (^^l{z)ip^p{w - ujf3)iiPa{z - w)(j){h,w) + ip^{z + w)(t){h, -w))- 
- u:^)^l{w){^^p{z - w)<i){h, w) + + w)<i){h, -u-))) = 
i[S^, Sp]+ (2rm-^^ai^HHV^»(' 4{Z)^1{W)^'1,{Z - W)) + 

2r^{h){-^l{z)^^p{w)^l{z + w) + ^%z)^l{w)^^p{z + w))) + 

Ai'a, i>l3]+ (20"'(^)(-9'a(^ - ^a)fpiw - UJp)ip^{z - w) + 

2(t)-'"{h){-^l{z - uJc.)v^p{w - ujp)^l{z + w)+ (^^(z - ujp)^l{w - W„)V9^(Z + w)) 



(D.l) 



where (a,/3, 7) is equivalent to (1,2,3) under cyclic permutations. The expression behind 
[i>„, vanishes due to ()('.12|) . At the same time the expression behind [v^, Sq\ has a pole at 
w = oj.y different from those of behind expressions [S^, So] and Sp\j^. Thus [i)^, Sq\ = 0. At 
the moment we have 

[S^, So] {2ip^{z)c^\w){<l)'^{h)^^{z -w) + ^~'"ih)^^iz + w)) - 
24>^{z)(p^{w){^'"{h)(t>''{z -w) + ^~'"{h)^^{z + w))) 

AiVaSp [^Pa{z - LOa)>f^p{w){ipl{z - w)(j){h, w) + Lp^{z + w)(j){h, -w))- 
(p^{z)ip^{w - UJa)i^'^iz - W)<i){h, W) + + W)(l}{h, -W))^ + 

AiDpS^ (ipl{z)^^^{w - ujp){ipl{z - w)<^{h, w) +ipl{z + w)<t){h, -w))- (^-2) 
^^p{z - up)^l{w){^^p{z - w)(t){h, w) + ip'^iz + w)(l){h, -w))^ = 
i[S^, Sp]+ {-2^l{z)^%w){r{ri)vl{z -w) + r'"{h)^^{z + w)) 
2^^p{z)^l{w){r{h)v%z -w) + r^{h)^%z + w))) 

Using then (|(].4|) and cancelling 2(j){z, K) we have 

[S^, So] {<f^{z)<i)f'{w){ip^{z -w) + if^iz + w))- 
ip^{w){<t)^{w)(t)^{z -w) + <j)^{-w)4>\z + w))) + 

^%z)^l{w - i^a){^p{z -w)+ ^p{z + w;))) + (0-3) 

2iDpSa (^^p%{z)ip^{w - UJp){ipaiz - w) + ipaiz + w))- 
ip^p{z - UJli)ip^{w){Lpii{z -w)+ iffsiz + w)) 
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= i[Sa, Sf3]+ (^-ip^{z)ip^{w){ipaiz - w) + Lfaiz + w))+ ip'^{z)cp^{w){ip(3{z - w) + ipp{z + 

To get the final result one should compare the structure of poles (tf = and z = —w). (|C.9IC.13|) . 
Other types of expressions can be simplified in the same way through the use of (|C.9IIC.13|1 . 
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